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Regimes of chemical reaction wave propagating in reactive gaseous mixtures, whose chemistry is 
governed by chain-branching kinetics, are studied depending on the characteristics of a transient 
thermal energy deposition localized in a finite volume of reactive gas. Different regimes of the 
reaction wave propagation are initiated depending on the amount of deposited thermal energy, 
power of the source and the size of the hot spot. The main parameters which define regimes of the 
combustion waves facihtated by the transient deposition of thermal energy are: acoustic timescale, 
duration of the energy deposition, ignition time scale and size of the hot spot. The interplay between 
these parameters specifies the role of gasdynamical processes, the formation and steepness of the 
temperature gradient and speed of the spontaneous wave. The obtained results show how ignition 
of one or another regime of combustion wave depends on the value of energy, rate of the energy 
deposition and size of the hot spot, which is important for the practical use and for risk assessment. 



PACS numbers: 47.70.Pq, 82.33.Vx, 47.40.Rs 
I. INTRODUCTION 

The initiation or ignition of a chemical reaction is one 
of the most important and fundamental problems in com- 
bustion physics. One needs to know how combustion 
starts and how the transient energy deposition influences 
the regime of the reaction wave which propagates out 
from a finite volume of reactive gas where a transient 
thermal energy were deposited - the hot spot. What 
type of combustion wave is formed depending on: the 
amount of energy actually added to a finite volume of 
reactive gas on a specific time scale, the power deposi- 
tion, and the ignition conditions, e.g. size of the hot 
spot, initial pressure, etc.? Long ago Oppenheim and 
Soloukhin^ recognized the importance of these concepts 
with the remark " Gasdynamics of Explosions is best de- 
fined as the science dealing with the interrelationship be- 
tween energy transfer occurring at a high rate in a com- 
pressible medium and the concomitant motion set up in 
this medium". The community of scholars has sought 
to address this perspective for many years represented 
by a vast combustion science literature too extensive to 
enumerate here. 

Transient thermal energy deposition into a reactive gas 
provides a source for ignition of either defiagration or 
detonation. Sufficiently fast and large energy addition 
can facilitate direct initiation of detonation. However 
the particular mechanism of the direct initiation of det- 
onation can be different. Detonation can be initiated 
by a strong shock (strong explosion) or it can arise as a 
result of the formation of an appropriate temperature 
gradient through the Zeldovich' gradient mechanism^. 
In most practical cases ignition arises from a small vol- 
ume of combustible mixture which is locally heated by 



energy input by means of an electric spark, hot wire, 
focused laser light and other related external sources. 
Such a transient energy addition process can generate 
a wide range of gas expansion processes depending on 
the amount and the rate of energy actually added and 
may result in the formation of the initial non-uniform 
distribution of temperature (see for example, Kassoy^ 
for the diverse range of fluid responses to localized, spa- 
tially distributed, thermal power addition into an inert 
gas). An example of an initial nonuniform distribution 
of temperature arises from the energy deposition of a 
spark-plug in an engine combustor— . In the general case 
it can be non-uniform distributions of temperature, pres- 
sure and/or concentration of reactants which determine 
further evolution of the reaction wave depending on the 
mixture reactivity and initial pressure. An example of 
concentration non-uniformity is a hydrogen gas leakage 
and its nonuniform distribution by convective mixing in 
a room. In all cases, a reaction wave arises from the in- 
duction time non-uniformity via the thermal explosion. 

The ignition problem is important for improving com- 
bustion safety and risk assessments of processes where 
hydrocarbons are oxidized at different initial conditions 
(concentration, temperature and pressure). How can we 
minimize "accidental" explosions in mines, chemical in- 
dustry and nuclear power plants? An important prob- 
lem of "hydrogen safety" is connected with leakage of 
hydrogen gas, subsequent mixing with air and the mix- 
ture explosion due to local heat release. It is worth not- 
ing that the problem in question is also of great interest 
for hydrogen storage, transportation and utilization and 
for the design and operation of perspective pulse detona- 
tion engines and homogeneous charge compression igni- 
tion (HCCI) engines. 
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For the first time possible regimes of chemical reaction 
wave ignited by the initial non-uniform distribution of 
temperature have been studied by Ya. B. Zeldovich using 
a one-step chemical reaction model^. The basic idea of 
the Zeldovich's concept was that a spontaneous reaction 
wave can propagate through a reactive material along a 
spatial gradient of temperature VT(x), with the velocity 

Usp = \{dnnd/dx)\-^ = \{dnnd/dTy\dT/dx)-^\ (1) 

where Tind{T(x)) is the induction time. The value of Usp 
depends only on the steepness of the temperature gra- 
dient. Then the regime of the formed combustion wave 
depends on the value of spontaneous wave velocity com- 
pared to the sound speed. 

Recently regimes of chemical reaction wave propaga- 
tion initiated by initial temperature non-uniformity in 
gaseous mixtures, whose chemistry is governed by chain- 
branching kinetics, were studied using a multispecies 
transport and detailed chemical model^. Possible regimes 
of the reaction wave propagation were identified for the 
stoichiometric hydrogen/oxygen and hydrogen/air mix- 
tures in a wide range of initial pressures and temperature 
level depending on the initial gradient steepness. 

The question that still remained unanswered is how 
the temperature gradient in the Zeldovich's concept of 
the spontaneous reaction wave^ arises. Kassoy and co- 
authorsi^— use a one step chemical reaction model to 
study how a temperature distribution adjacent to a pla- 
nar boundary, generated either by direct deposition of 
transient, spatially distributed thermal power into a de- 
fined volume of reactive gas, or by conduction through 
the boundary into the gas, leads to planar detonation 
initiation. The authors seek to understand the magni- 
tude of the energy addition deposited on specific time 
and length scales required to produce conditions that will 
lead to detonation initiation. The important difference 
between a one-step and detailed chemical models is seen, 
for example, from the result of^, where it was shown the 
ignition energy for methane-air computed using a one- 
step model differs by two orders of magnitude from the 
experimentally measured value. 

Liberman et. al^ take a different approach to show that 
steepness of an imposed temperature gradient (the length 
scale at fixed temperature difference) required for initi- 
ating combustion regimes, in particular a detonation, in 
gaseous mixtures with chain-branching kinetics, may dif- 
fer up to several orders of magnitude from that obtained 
using a one step global chemical reaction model. The en- 
ergy of ignition for hydrogen-oxygen gaseous mixture has 
been calculated in^^ using a detailed reaction mechanism 
and a multispecies transport model, but difference in ig- 
nition of different combustion regimes remained unan- 
swered. 

The purpose of the present paper is to study how the 
initial temperature gradient is formed depending on the 
rate and amount of energy addition, on the size of hot 
spot and on the gasdynamic processes emerging in the 
region of transient energy deposition. The problem of 



great practical importance is: how is temperature non- 
uniformity initiating different combustion regimes arise? 
How do the gasdynamic processes caused by the energy 
addition influence the formation of the temperature gra- 
dient and how it depends on the rate and amount of en- 
ergy input and on the size of hot spot? Solution of this 
problem answers the question of great practical impor- 
tance: what is the amount of energy and how it should 
be deposited to ignite certain regimes of combustion 
wave. The paper presents new results on classification 
of the propagation regimes of chemical reaction wave ini- 
tiated by the transient energy deposition in gaseous mix- 
tures using high resolution numerical simulations of re- 
active Navier-Stokes equations, including a multispecies 
transport model, the effects of viscosity, thermal conduc- 
tion, molecular diffusion and a detailed chemical reac- 
tion mechanism for hydrogen-oxygen mixture which is 
the quintessential example of chain branching reactions 
whose chemical kinetics is well understood and whose 
detailed chemical kinetic models are well known and rel- 
atively simple. Such a level of modeling allows clear un- 
derstanding of the feedback between gasdynamics and 
chemistry, the principal point when studying unsteady 
process of ignition, not easily be captured using simpli- 
fied gas-dynamical and chemical models. 

High fidelity reliable numerical simulations of the 
present study are performed to identify ignition processes 
in homogeneous hydrogen-oxygen mixtures caused by the 
localized energy deposition. The distinct ignition regimes 
are identified depending on the size of the hot spot, en- 
ergy amount, the duration of energy addition and initial 
pressure. The obtained results open perspectives for un- 
derstanding of how to avoid or on the contrary what are 
the conditions to initiate different combustion regimes 
(slow deflagration, fast deflagration, detonation). The 
performed analysis reveals the conditions when detona- 
tion is initiated as a result of direct initiation by a strong 
shock wave, or when it results from the formation of suit- 
able temperature gradient for the Zeldovich's mechanism 
of detonation triggering. The amount and the rate of 
energy addition needed to form a proper gasdynamics 
nonuniformity for triggering either detonation or defla- 
gration waves are found. 



II. PROBLEM SETUP 

We consider uniform initial conditions and a transient 
external source of energy localized on the scale of the 
"hot spot" < a; < L, where energy Qig is added during 
the time Atg . Gasdynamics of the explosion initiated by 
the localized energy deposition is characterized by the in- 
terrelationship between the rate of energy transfer with 
the time characterizing energy deposition Atg in the hot 
spot of size L and the characteristic times of the problem. 
For the sake of simplicity we assume that the rate of the 
energy addition is constant in time, so that total energy 
deposition into the hot spot is Qig = WAtg, where W 
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is the power of the external source of energy, and Aig 
is time of the hot spot transient heating. The character- 
istic acoustic time ta = L/a{T) defines the concomitant 
motion setup in the gaseous mixture, where a{T) is the 
speed of sound. If Atq ^ ta local heat addition occurs 
as a nearly constant volume process and the tempera- 
ture elevation within the hot spot is accompanied by a 
concomitant pressure rise. Subsequent expansion of the 
hot spot driven by the large pressure gradient between 
the hot spot and the ambient gas causes compression 
and shock waves in the surrounding ambient gas. When 
Atq ^ ta the acoustic waves have enough time for pres- 
sure equalization and the local heat addition occurs at 
nearly constant pressure. 

The local heat of gaseous mixture within the hot spot 
leads to the heat being concentrated in a small region 
from which heat propagates into the surrounding gas ac- 
cording to well known solution of the equation of thermal 
conduction^®. For the one-dimensional case the temper- 
ature distribution (for constant thermal diffusivity) is 

T(x, t) (X exp(~xV4xi) (2) 

where x = i^{T)/ pCp is the coefficient of thermal diffu- 
sivity. It follows from Eq. ^ that heat propagates at the 
distance x ^ \/4xi and temperature in the surrounding 
gas increases noticeably due to propagation of a thermal 
wave during the time 

tT ^ x^/x (3) 

At the very beginning the thermal wave propagates 
with the velocity dx/dt ~ \/xA overtaking the shock 
wave, which propagates with a velocity approximately 
equal to the sound speed in the heated gas. 

However, very soon the shock wave overtakes the ther- 
mal wave, so that characteristic thermal wave time is 
effectively much longer than the acoustic time. For ex- 
ample, tT ~ Ix ~ 0-1^ and ta ^ 2/iS for x ~ 1mm for 
hydrogen- oxygen mixture at P = latm. 

Whether a chemical reaction starts and which pro- 
cesses define the regime of combustion wave depend on 
the interrelationship between Atg, ta, tT and induction 
time tind{T) at the existing temperature and pressure. 
The induction time is the time scale for the stage of en- 
dothermic chain initiation and branching reactions (in 
the case of a global one-step reaction this is the time- 
scale for the maximum reaction rate). Here it is suit- 
able to use the scale of ignition time tign , characterizing 
the length of induction phase after or during the tran- 
sient energy deposition tind{T, P). The induction time is 
measured experimentally and determines local properties 
of the combustible mixture depending on its thermody- 
namic state. Dependence of the induction time on tem- 
perature for a hydrogen- oxygen mixture at initial pres- 
sures Po = latm and Po — lOatm is shown in Fig. [T] 

For the temperature range T = (1100— 1500)K, where 
the exothermic reaction starts, the induction time is 




T(K) 

FIG. 1. Induction time for hydrogen-oxygen stoichiometric 
mixture at different temperatures and pressures: Po = latm 
(solid line) and Po = latm (dashed line). 



about ten microseconds. If the reaction has started, fur- 
ther heating and energy deposition do not matter and do 
not influence the formed combustion regime. Thus, pri- 
mary interest is focused on the regimes of ignition with 
Atq < tign- For a very short time of energy addition, 
much less than the acoustic time, the mixture in the hot 
spot can be heated to any temperature and the ignition 
regime will be determined by the induction time at that 
temperature and accompanying pressure. In case of a 
more extended energy deposition, the ignition regime will 
depend on the size of the hot spot (and correspondingly 
on ta) and the relation between Atq and ta as mentioned 
earlier in reference to the features of the hot spot expan- 
sion process. 

For a one-dimensional formulation the hot spot repre- 
sents a region < x < L, where energy is deposited with 
the rate W{t) — dQ{t)/dt. The specific internal energy of 
the mixture is changed as Q{t) = Q{t — At) + AQ within 
the hot spot during time interval At. Total energy at the 
end of the energy addition depends on the gas-dynamic 
motion 



Atq L 

Qtota; = ^"3 • /^(^' • ^2;. (4) 

t=0 x=0 

Total energy transfer till the formation of a 1-D sta- 
tionary combustion wave far from the hot spot (at x 3> 
L) can be viewed as the ignition energy required for ini- 
tiation of the particular combustion regime. 

The governing equations are the one-dimensional time- 
dependent, multispecies reactive Navier-Stokes equations 
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including the effects of compressibility, molecular diffu- 
sion, thermal conduction, viscosity and chemical kinetics 
with subsequent chain branching, production of radicals 
and energy release. 
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The initial conditions at i = are constant pressure 
and zero velocity of the unburned mixture. At the left 
boundary at x = the conditions are for a solid reflecting 
wall, where u(0, t) — and the initial temperature T{t = 
0) = To. 

Here we use the standard notations: P, p , u, are pres- 
sure, mass density, and flow velocity, Yi ~ Pi/p - the 
mass fractions of the species, E = e ^ /2 - the to- 
tal energy density, e - the internal energy density, Rb 
- is the universal gas constant, rrii- the molar mass of 
i-species, Ri = RB/nii, n - the molar density, Uij- the 
viscous stress tensor, c„ ~ "^CviYi - is the constant vol- 



ume specific heat, Cyi- the constant volume specific heat 
of i-species, hi - the enthalpy of formation of i-species, 
K,{T) and /i(T) are the coefficients of thermal conductiv- 
ity and viscosity, Di(T) - is the diffusion coefficients of 
i-species, {dYi/dt)^f^ - is the variation of i-species concen- 
tration (mass fraction) in chemical reactions. 

The equations of state for the reactive mixture and 
for the combustion products were taken with the tem- 
perature dependence of the specific heats and enthalpies 
of each species borrowed from the JANAF tables and 
interpolated by the fifth-order polynomial o^^i^^ . The vis- 
cosity and thermal conductivity coefficients of the mix- 
ture were calculated from the gas kinetic theory using the 



Lennard- Jones potential^^ Coefficient of the heat conduc- 
tion K = fiCp/ Pr for the mixture as a whole is expressed 
via the viscosity p and the Prandtl number, Pr — 0.75. 

The numerical method is based on splitting of the Eu- 
lerian and Lagrangian stages, known as coarse particle 
method (CPM)2S. A detailed description of the modified 
CPM optimal approximation scheme, details of the equa- 
tions and transport coefficients and the reaction kinetics 
scheme together with the reaction rates appear in^L2^. 
The numerical method is thoroughly tested and success- 
fully used in many practical applications^^"^'*. 

The convergence of the solutions is of paramount im- 
portance to verify that the observed phenomena are suf- 
ficiently well resolved, especially when CFD simulations 
are used with a detailed chemical reaction model. The 
convergence and resolution tests have shown that the 
resolution of 50 computational grid cells over the width 
of a laminar flame (for example, with the grid cell size 
A = 0.0064mm at Pq = latm, when the width of a 
laminar front is 0.24mm, and much smaller for higher 
pressure) provides sufficiently good convergence and cor- 
rectly captures the details of the observed processes (see 
Appendix in Refi^). 



III. COMBUSTION REGIMES: RAPID ENERGY 
DEPOSITION - MICROSECOND TIME SCALE 

We consider combustion regimes initiated by energy 
deposition at the initial pressure and temperature Pq = 
latm and Tq = 300K in a hydrogen-oxygen stoichiomet- 
ric mixture. First we consider cases when the time scale 
of the energy deposition is comparable to or shorter than 
the acoustic time scale and less than the induction time at 
the ignition temperatures. The acoustic time is ta ~ 20/iS 
for the hot spot of size L — 1cm, and ta ~ 2/is for 
L = 1mm (ao(r = 300A') = 539m/s). Rapid energy ad- 
dition into the hot spot on time scales much shorter then 
acoustic time causes almost uniformly fast elevation of 
pressure and temperature resulting in the volume explo- 
sion. Figure [2] shows the calculated temporal evolution 
of temperature at the center of the hot spot {x — 0) for 
different values of the transmitted energy, when the en- 
ergy addition time is very short: Atq O.lps < ttgn ^ ta- 
Temperature and pressure of the mixture in the hot spot 
depend on the energy transmitted to the hot spot. After 
the end of the energy deposition process, the induction 
period reaction starts. After about lOps the stationary 
combustion regimes are established (see Fig. [5]). It should 
be noted that for each value of deposited energy there is 
some definite temperature and pressure at which the re- 
action starts and the combustion regime is produced by 
the volumetric explosion at these conditions. 

For a less rapid process of energy deposition {Atq = 
5ps < tign < ta), a large pressure jump is formed at the 
boundary of the hot spot. If the power is large enough, 
subsequent events represent the decay of the initial dis- 
continuity consisting of a compression wave propagating 
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FIG. 2. Temperature evolution in the center of the hot spot 
for rapid submicrosecond energy deposition Atg = (0.1 
0.2)/is; 1 - deflagration {Q — 1.9kJ /m^), 2 - deflagration 
(Q = 2AkJ/m^ ), 3 - detonation (Q = i.QkJ/m^ ). 




FIG. 3. Evolution of temperature profiles illustrating det- 
onation formation in the energy release region. L = 1cm, 
Atg = 5/is. Profiles are presented for time instants with in- 
terval Atg = 0.5/iS. 




FIG. 4. Evolution of pressure profiles illustrating detonation 
formation in the energy release region. L = 1cm, Atg = 5/is. 
Profiles are presented for time instants with interval Atg = 
0.5/is. 



in the direction a; > 0, which steepens into the shock 
wave and the rarefaction wave propagating to the left 
from the boundary of hot spot (in the direction x — Q) 
with the velocity equal to the local sound speed. Such a 
scenario, which is similar to a strong point explosio n^^i^^ 
results in the direct triggering of a detonation wave if a 
concomitant shock wave at the right boundary is strong 
enough. A more interesting scenario emerges in the case 
of a weaker shock and is shown in Fig. [3] and Fig. |4l 
These figures show the calculated transient evolution of 
the temperature (Fig. [3]) and pressure (Fig. 2]) profiles 
inside the hot spot (L = 1cm) during the energy deposi- 
tion. The rarefaction wave propagating to the left creates 
shallow temperature and pressure gradients on the scale 
of about the size of the hot spot. At an initial pressure 
of Pq = latm the temperature gradient with the temper- 
ature difference and length scale L ^ 1cm as in Figs.|3l|4] 
cannot trigger detonation. However, since the pressure 
of the heated mixture increased during the heating up to 
P « 4atm, this temperature gradient can produce a deto- 
nation through the Zeldovich gradient mechanism. Time 
evolution of the temperature and pressure profiles shown 
in Figs. [21 m demonstrate also the emergence of the spon- 
taneous reaction wave and its coupling with the pressure 
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FIG. 5. Temperature evolution in the center of the hot spot 
for case of low energy deposition {ta < Atq < tind)- L — 
1mm, Atq — 16/is (solid line), Po — lOatm (dashed line). 
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FIG. 6. Evolution of temperature profiles in the hot spot 
during energy deposition (solid lines) and during induction 
phase (dashed lines). L = 1cm, Atg = 16/is. Profiles are 
presented for time instants with interval At = l/is . 



wave leading to the detonation initiation through the Zel- 
dovich mechanism. 

In case when the ignition time is greater than acoustic 
time and the energy addition time is less than acous- 
tic time Atq < ta < ttgn the gradient induced by the 
rarefaction wave forms on the stage after the end of en- 
ergy addition. During the relatively long induction phase 
the acoustic perturbations equalize pressure in the en- 
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FIG. 7. Evolution of temperature (dashed lines) and pressure 
(solid lines) profiles illustrating ignition process in the energy 
release region, L = 1cm, Atg = 16/is. Profiles are presented 
for time instants with interval At — Ifis . 



ergy deposition zone and further ignition processes evolve 
at constant pressure from the steady temperature gradi- 
ent. Temperature in the top of gradient remains nearly 
constant till the ignition takes place (see Fig. [5]). The 
combustion regime forming in this case depends on the 
environmental conditions (e.g. pressurSi^): at Pq = latm 
such a temperature gradient causes deflagration wave for- 
mation, at Po = lOatm it causes detonation wave. 

An interesting scenario takes place when the energy 
deposition time is slightly longer then acoustic time, 
Atq = 20/is, so that ta < Atq < tign- At the beginning 
the gasdynamics process is similar to the previous case. 
However, the heating time is longer, and the rarefaction 
wave has time to reach the left wall, to be reflected from 
the wall and return to the edge of the hot spot during 
time of the energy deposition (15 in Ref. — ). As a result 
even before the reaction started the hot spot expands 
with the temperature profile in the mixture consisting of 
a shallow temperature gradient in the direction to a; > 
and a steep temperature gradient in the direction X = 0. 
The initial stage of the temperature time evolution, be- 
fore the reaction has started, is depicted in Fig. [51 By the 
time the reactions begins at the top of the gradients, the 
right side gradient is too steep to trigger a detonation. 
Instead a deflagration wave propagating to the right is 
ignited Ref. — ). At the same time due to the elevated 
pressure the left shallow gradient can facilitate trigger- 
ing detonation through the Zeldovich gradient mecha- 
nism. The corresponding time evolution of temperature 
and pressure profiles showing propagation of the sponta- 
neous wave to the left along the temperature gradient, its 
coupling with the pressure waves and detonation initia- 



tion as well as deflagration initiation by the right steep 
temperature gradient are depicted in Fig. [T] 

The rapid energy deposition such that the heating 
time, Aig, is comparable to the characteristic acoustic 
time scale of the volume, ta^ always results in the shock 
waves propagating away from the hot spot. A particu- 
lar scenario of the resulting combustion regime depends 
on the size of the hot spot, though the basic physics ap- 
pears to be similar to that described above. In the case 
of the smaller size of the hot spot (L = 1mm) scenario 
of the combustion regime ignition may differ because the 
scale of the temperature gradient created by the rarefac- 
tion wave may not be compatible with the detonation 
formation in real mixtures. There are many different 
scenarios that include direct detonation formation by a 
strong enough shock wave in the context of a thermal 
explosion, or the shock waves propagating away from the 
hot spot producing ignition of the fast deflagration prop- 
agating behind the shock waves (regime 3 according to 
the classification in Refj^). 



IV. COMBUSTION REGIMES: MILLISECOND 
TIME SCALE OF ENERGY DEPOSITION 

Sufficiently rapid and large amount of thermal energy 
deposition into a reactive gas can trigger either direct ini- 
tiation of detonation through a constant volume explo- 
sion, or through the Zeldovich gradient mechanism due 
to the shallow gradient formed by the rarefaction wave 
at the increased pressures in the hot spot region. In case 
of the rapid but relatively small energy deposition the re- 
sulting regime is a fast deflagration wave Ref<^. The sce- 
nario for low power thermal energy addition over a longer 
period of time is different. If the acoustic time is much 
less than the energy deposition time, ta <S^ Atg < tign 
, then there is enough time for pressure to be spatially 
homogenized by acoustic waves. In this case there are no 
strong compression waves emitted from the hot spot, and 
the combustion regimes initiated by the energy deposi- 
tion into the hot spot depend essentially on the steep- 
ness of the temperature gradient, which is formed by the 
thermal wave and gas expansion in the vicinity of the hot 
spot. 

During the time of energy deposition Atq the thermal 
wave propagates away from the hot spot at the distance 
XT/mm = (xAtQ)i/2 « 0.9{AtQ/msy^^ . Some of the 
mass in the volume of the hot spot heated by the added 
energy, therefore of higher temperature, flows away as the 
temperature increases and the density falls. The expelled 
mass together with the thermal wave give rise to the tem- 
perature gradient in the surrounding mixture behind the 
boundary of the hot spot. The temperature profile is al- 
most linear because of the weak temperature dependence 
of the coefficient of thermal conduction {k oc T°'^^). 
Fig. [5^ shows the temperature gradient formed at the 
end of the energy deposition for Atq = 1000/xs, L = 1mm 
(solid lines) . A greater distance compared to that created 
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FIG. 8. (a) Evolution of temperature profiles in the hot spot 
during energy deposition (solid lines) and after it during the 
combustion wave formation (dashed lines). (b)Evolution of 
H2O concentration (dashed lines) and ff-radical concentra- 
tion (solid lines) profiles illustrating combustion wave forma- 
tion on the gradient formed in the energy release region.!/ = 
1mm, Atq = 1000/is. Profiles are presented for time instants 
with interval At — 5fi s. 



by the thermal wave alone is due to the hot spot expan- 
sion, the decreased density and increased temperature 
in the hot spot during the process of energy deposition. 
However, according to Ref.^ this type of linear tempera- 
ture gradient is too steep to initiate detonation through 
the Zeldovich gradient mechanism and as a result a de- 
flagration wave is initiated. The thermal wave and the 
gas expansion are too slow to expand temperature and to 
form a temperature gradient compatible with the detona- 
tion formation in real mixtures at atmospheric or lower 
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FIG. 9. Temperature evolution in the center of the hot spot 
for case of slow energy deposition {ta <C Atg). L = 1mm, 
Po = latm (solid line), Po = lOatm (dashed line). Atg = 
1000/is for Po = latm and 1100/is for Po — lOatm. 



pressures Ref.^. Long before the thermal wave moves 
away a sufficiently long distance the temperature of the 
mixture rises to ignite the reaction, so that either a de- 
flagration wave or a fast deflagration wave are initiated 
according to the classification of Refi^. In Fig. [8^ dashed 
lines show the deflagration wave formation out from the 
formed temperature gradient. Fig.lSb represents the evo- 
lution of i?-radical and H2O concentration profiles while 
the combustion wave is forming. 

While speed of sound does not depend on pressure, 
the induction time tind{T) at the temperature range 
(1100-^ 1200)K is considerably longer at larger pressures 
(see Fig. [ij. This leads to a significant delay of ignition 
for the same energy deposition regime at higher pres- 
sure. As an example, the growth of temperature at the 
center of the hot spot L = 1mm and at instants when 
the reaction starts are shown in Fig. 9 for initial pres- 
sures Pq = latm and Pq = lOatm and for Atg = 1ms. 




2 

X (mm) 

FIG. 10. (a) Evolution of temperature (dashed lines) and 
pressure (solid lines) profiles illustrating detonation forma- 
tion on the gradient formed in the energy release region, (b) 
Evolution of H2O concentration (dashed lines) and //-radical 
concentration (solid lines) profiles illustrating detonation for- 
mation on the gradient formed in the energy release region. 
L — 1mm, Atg — 1100/is, Po = lOatm. Profiles are presented 
for time instants with interval At = 2/xs. 



Since the coefficient of thermal conductivity does not 
depend on pressure, and the steepness of the temperature 
gradient for direct detonation initiation through the Zel- 
dovich gradient mechanism decreases considerably with 
the increase in pressure Ref.^, the temperature gradient 
created by the thermal wave can trigger detonation at 
high enough initial pressure. To elucidate the process we 
consider a relatively small hot spot of size, L = 1mm 
at initial pressure Pq = lOatm. For a sufficiently long 
energy deposition time, the thermal wave creates gra- 
dient appropriate for the detonation initiation through 
the Zeldovich mechanism at Pq = lOatm Ref.^ (dashed 
line in Fig. [5]). Figure [TUk shows the calculated temporal 



evolution of temperature and pressure profiles illustrat- 
ing formation of the temperature gradient outside of the 
hot spot, L = 1mm, the development of the spontaneous 
wave along the gradient and transition to detonation for 
the energy deposition time Atq = 1ms. Fig. [TUb shows 
the corresponding evolution of the concentration profiles 
for //-radicals and the combustion products {H2O). 
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FIG. 11. Temperature evolution in the center of the hot 
spot for case of rapid energy deposition {ta > Atg) in one- 
dimensional case (dashed line) and three-dimensional case 
(solid line), L = 1mm. 



FIG. 12. Temperature evolution in the center of the hot 
spot for case of rapid energy deposition {ta > Atg) in one- 
dimensional case (dashed line) and three-dimensional case 
(solid line); L — 1mm. 



V. ENERGY OF IGNITION 

The results obtained clearly allow us to estimate the 
energy required for ignition of a particular combustion 
regime. At the same time it should be noted that the 
amount of the ignition energy obtained by extrapolating 
results of the one-dimensional problem most likely will 
not match the actual value of the ignition energy for the 
three-dimensional problem where the process is associ- 
ated with a three-dimensional expansion and converging 
rarefaction wave. This difference is particularly impor- 
tant for the initiation of detonation. In this case the three 
dimensional expansion additionally enhances the rarefac- 
tion, leading to less suitable conditions for a detonation 
initiation. A large number of 3D simulations of ignition 
due to the energy addition to the spherical hot-spot were 
used to verify realization of different combustion regimes, 
and to compare the ignition energy obtained from 3D 
spherically symmetric model to that extrapolated using 
the ID model and to assess effect of spherical expansion 
on the ignition process. 

Figure [TT] shows the temperature evolution in the cen- 
ter of the hot spot for rapid energy deposition (ta ^ Atg) 
in a one-dimensional case where L ~ 1mm and in the 
three-dimensional spherical hot spot of radius i? — 1mm. 
In both cases the final temperature in the center of the 
hot spot tends to the same value. Temperature oscil- 
lations in the spherical hot spot are caused by conver- 
gence of the rarefaction wave to the center and reflec- 
tion, which is enhanced by the expansion of the gaseous 
spherical hot spot, where the density decreases approxi- 



mately as poc ^1 — (r/i?)^^ . Since the rarefaction is 

much stronger for the 3D spherical expansion, the final 
temperature of the spherical target is considerably lower 
than in the planar case at the same parameters of energy 
deposition if the rarefaction refiects before the energy is 
deposited it a < Atg). For the same reason, the actual 
energy required to initiate a detonation is larger for the 
spherical target than that in the planar case. Even if the 
temperature of the hot spot rises high enough for initi- 
ating detonation in the planar case, it is reduced by the 
converging and refiecting rarefaction wave and detona- 
tion can not be ignited. An example of such a scenario 
is shown in Fig. [T^l where for the same conditions the 
detonation is initiated in the plane geometry and it is 
not for a spherical target. 

The larger the acoustic time ta compared to the time 
of energy deposition the less the influence of the rarefac- 
tion wave on the detonation initiation. This means, for 
example, that with the increase of the hot spot size (note, 
that ta oc L) the initiation of detonation requires less en- 
ergy deposition into the specific volume of the hot spot 
for a given time of power deposition Aig, although the 
total deposited energy will be larger than in the case of 
detonation initiation in a smaller hot spot using a higher 
level of power deposition. This tendency is illustrated 
in Fig. 1131 For a shorter period of power deposition the 
energy amount per unit of volume capable for detona- 
tion initiation decreases. The lower limit for hydrogen 
detonation initiation can be obtained using a short sub- 
microsecond laser pulse focused in a sub-millimeter area 
(see e.g. Rei^) and it is estimated as ~ 10~^mJ, which 
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FIG. 13. The minimal energy required for detonation initia- 
tion in H2-02 depending on the size (radius) of the hot spot 
for Atq = 1/is. SoUd Une is total deposited energy; dashed 
line is the specific energy {Q/V[m3 ■ m^). 
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transient energy deposition into a stoichiometric mixture 
of hydrogen-oxygen leading to the ignition of different 
regimes of combustion. It is shown that depending on 
the parameters of energy deposition (deposited energy 
amount, deposition time scale and size of the hot spot) 
there are two main mechanisms of reaction wave initi- 
ation: the Zeldovich gradient mechanism Ref.^ and the 
volumetric thermal explosion (which actually represents 
one of the asymptotics of the Zeldovich mechanism for 
the gradient of zero steepness) . For practically important 
time scales the principal scenarios of ignition are: 1) for 
sub-microsecond pulses the volumetrical explosion takes 
place inside the hot spot; 2) for microsecond pulses the 
gradient of temperature and pressure arises on the pro- 
file created by the rarefaction wave and ignition starts 
via Zeldovich mechanism on the gradient of induction 
time; 3) for millisecond pulses gasdynamical expansion 
gives rise to a temperature gradient at approximately 
constant pressure, and the ignition starts according to 
the Zeldovich mechanism Refi^ with all the features in- 
herent to chain-branching chemistry disclosed in Ref^^. 
In the three-dimensional case spherical expansion of the 
hot spot weakens the generated shock wave in favor of an 
intensified rarefaction wave. It results in sufficient drop 
in temperature and pressure in the hot spot on the time 
scales of the order of acoustic time. Thus, for the same 
conditions as in one-dimensional less intensive 

combustion regime arises. The deflagration regimes are 
less sensitive. However to obtain detonation one should 
increase sufficiently the power of the energy source. For 
example for L — 1mm and Afg = l//s the energy amount 
is about 10-12 times larger compared to that obtained via 
extrapolating of the results of ID model. This is clearly 
seen from Fig. 1141 which shows the ratio of energy re- 
quired for the detonation initiation in ID and 3D cases. 
The calculations were done for the fixed time of the en- 
ergy deposition Atg = 1/xs, for the planar and spherical 
hot spots of different sizes ( i in ID case and i? -f L for 
3D case). In particular, it is seen that with the increase 
of the hot spot size, and corresponding increase of the 
acoustic time i^, the role of a rarefaction wave becomes 
less important. 



FIG. 14. The ratio of energies required for detonation initi- 
ation in H2 — O2 obtained for 3D and ID model versus the ACKNOWLEDGMENTS 
size (radius) of the hot spot for Afg = 1/i. 



agrees with the extrapolation of the dependence shown 
in Fig. [131 

VI. DISCUSSION AND CONCLUSIONS 

In the present paper we used detailed chemical kinetics 
and transport models to study consequences of localized 
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